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Abstract 

A model of a high subsonic jet with a nearby array of flexible, aircraft-type panels is studied 
numerically in two dimensions. The jet is excited by a limited duration, spatially localized 
starter pulse in the potential core. The long time evolution of unsteady disturbances in 
the jet, the responses of the panels and the ensuing radiation are computed. The results 
show that the spectral response of both the jet and the panels is concentrated in a relatively 
narrow frequency band centered at a Strouhal number (based on jet exit velocity) of approx- 
imately 0.25 and associated harmonics. The loading on the panels generally increases with 
downstream distance. Panel radiation is weakest in upstream directions. Interior zones of 
silence, due to destructive interference of radiation from the panels, are observed. 
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1 Introduction 


This paper describes the results of a numerical simulation of jet noise in the presence 
of four flexible aircraft-type panels in a panel-stringer assembly. The simulation is based on 
a model which fully couples the fluid dynamics of the jet flow to the panel motion and the 
resulting acoustic radiation. The primary objective is to determine the role played by near 
sonic jet exit conditions on installation effects from the nearby flexible structure and on the 
response and the acoustic radiation from the structure. 

In this paper we consider the long time evolution of unsteady disturbances in a nearly 
sonic jet exiting from a converging nozzle together with panel responses and ensuing radi- 
ation from the panels. In previous work, jet acoustics, panel response and radiation have 
been considered for lower speed jets exiting from straight pipes, both at rest [1] and in for- 
ward motion [2]. The present work focuses on jet evolution over much longer timescales, 
emphasizing in particular vortex shedding from the nozzle lip. 

In previous analyses, the exact sources of jet noise have been identified from the basic 
equations of fluid dynamics [3, 4, 5, 6]. Generally, these exact sources have to be modeled 
in some way for computation. An important feature of the present method is the direct 
computation of at least some of the natural sources of jet noise, namely fluid dynamical 
instability waves which develop due to instability of the jet shear layer. Experiments have 
demonstrated the existence of large scale structures or instability waves in jets [7, 8, 9]. 
These structures are believed to act as sources of sound, a point also confirmed by analytical 
studies [10, 11, 12, 13] and computations [9, 14, 15]. 

The geometry of the computational model is shown in figure 1. The Euler equations are 
solved in two domains; the jet domain and the radiation domain separated by an array of four 
flexible panels. Panel response and radiation are also computed and are fully coupled to the 
fluid dynamics in the sense that at each timestep the fluid dynamics (Euler) computation 
provides the pressure difference across the panels, thereby allowing computation of panel 
displacement and velocity. The resulting panel velocity then serves as a boundary condition 
for the Euler computation. Thus no modeling of the panel excitation is needed. Other 
approaches to the general problem of describing panel loadings are possible, e.g., [16]. 
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The paper is organized as follows. In the next section there is a description of the 
model and a discussion of the numerical method and boundary conditions. Following 
that the results are presented first for the case of a free jet and then for the full jet 
flow/acoustic/structural interaction. Finally, the results are summarized and conclusions 
are drawn. 


2 Numerical Method 


The computational domain is shown in figure 1. Unsteady pressure, density and velocity 
are computed in two regions. The lower domain in figure 1 contains the jet, exiting from a 
contracting nozzle of width D, while in the upper domain, (modeling the aircraft interior) 
only disturbances due to the panel radiation are present. These two regions will be referred 
to as the jet and radiation domains, respectively. The wall boundary is a rigid wall con- 
taining four adjacent flexible panels (denoted as panels 1-4 in figure 1) with rigidly clamped 
boundaries. The panels vibrate in response to excitation by sound from the jet, thereby 
radiating acoustic disturbances into both domains. In the jet domain, the panel radiation 
level is significantly smaller than the acoustic disturbances generated in the jet, in contrast 
to the radiation domain where panel radiation is the sole source of acoustic disturbances. 

The numerical method involves coupling the computation of a nonlinear equation gov- 
erning the panel responses (the beam equation) to an Euler computation performed in both 
the jet and radiation domains. The panel vibration is fully coupled to the fluid dynamics in 
that at each time step the pressure difference across the panels, computed from the Euler 
computations, serves as a forcing term for the beam equation. Similarly, the displacement 
obtained from the beam equation is differentiated in time and is imposed as a boundary con- 
dition for the Euler computations. The numerical method has been described in detail [1]. 
Thus, the presentation here will be brief. 

The nonlinear beam equation is 


n d * 2 \t / \^ z , d 2 Z dz , 

Db d. ^ + pbh W + 7 dt ~ p 


(i) 


where 2 represents the beam transverse deflection, pb is the mass per unit volume of the 
beam, h is the beam thickness, 7 is the physical damping, Db = Mh 3 / 12(1 — v 2 ) is the 
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stiffness of the beam, M is the modulus of elasticity and v is the Poisson ratio of the beam 
material. The coefficient N x (z) of the nonlinear term represents the tension created by the 
stretching of the plate due to bending. The pressures in the radiation and jet domains are 
and respectively. The solution of equation (1) is obtained at each time step using an 
implicit finite difference method. The panels are assumed clamped at both ends. 

The coupling of the beam computation to the Euler computation occurs through the 
forcing term given on the right hand side of equation (1). The pressures p + and p~ are 
obtained from the Euler computation using an explicit scheme. The displacement at the new 
time level is then obtained from solving equation (1) one time step. The normal velocity, v, 
is then obtained from differentiating z and employed as a boundary condition to complete 
the update to the Euler computation. Since this procedure is employed at each time step, 
the fluid and structural calculations are fully coupled. 

The Euler equations are solved in conservation form for the vector 

w = (p, pu, pv, E) t , 

where p is the density, u. v are the x and y components of the velocity respectively and E is 
the total energy per unit volume, 

E — ±/>(u 2 + v 2 ) + CypT , 

where T is the temperature and is the specific heat per unit volume. The pressure, p , is 
obtained from the equation of state. The Euler equations are solved separately in both the 
jet and radiation domains. 

In the jet domain the Euler equations are modified to account for the jet flow. The jet 
exits from a nozzle of width D and the solution is computed both within and exterior to 
the nozzle. The Euler equations are modified to account for two different non-homogeneous 
forcing terms [1]. One term serves as a starter pulse to excite the jet. It corresponds to a 
localized source of mass injection at the location (x s , yj), where yj is the location of the jet 
axis (3.5 D from the wall) and x s is approximately 1.15 D. An alternative approach, involving 
time harmonic excitation of the jet is described in [17]. The second forcing term is designed 
so that in the absence of the starter pulse the solution to the Euler equations would be 
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a stationary profile corresponding to a spreading jet. The inclusion of this term separates 
the computation of the disturbance, in particular the resulting instability waves, from the 
computation of the mean flow (i.e. the spreading jet). Thus, the resulting system of equations 
allows for the simulation of instability waves and the resulting sound generation, together 
with the curving of acoustic waves in the jet flow field, without requiring the computation of 
the spreading jet itself. Although this is a simplified model, it captures many of the observed 
features of instability wave generated jet sound and permits high resolution computation of 
the coupling of jet noise with the flexible panels and the resulting radiation from the panels. 
In particular, the model allows for computation of the natural sources of jet noise (the 
instability waves) together with the sound radiated by these sources. 

The initial conditions are taken to be ambient data in the radiation domain and the mean 
state in the jet domain. The model for the mean axial and normal velocities is obtained 
from [18] in Cartesian coordinates. In contrast to previous calculations [1, 2] at lower jet 
speeds, the current study requires a mean flow model which addresses the compressibility of 
the flow. Thus, the mean pressure and density are also modeled. Inside the nozzle, the quasi 
1-D equations of isentropic flow are employed for pressure and temperature. Outside the 
nozzle, the jet pressure is obtained from [19]. Centerline temperature is assumed isentropic 
and off-centerline temperature variations follow the axial velocity relaxation. Density follows 
from the state equation both inside and outside the nozzle. It is imperative that the mean 
flow be smooth up to and including its second derivative for the numerical simulation of 
instability flows [20, 21]. Therefore, great care was taken in modeling the mean flow. 

The boundary conditions are as follows (refer to figure 1): 

1. Bounding wall - rigid conditions are imposed except for the flexible panels which 
are treated as described above. 

2. Nozzle exterior - impedance boundary conditions are used on the exterior of the 
nozzle, simulating the use of an absorbing material to absorb waves incident on 
the nozzle from the exterior. 

3. Nozzle interior - rigid conditions are used for the terminal two thirds of the nozzle. 
The first third (near the inflow boundary) employs impedance conditions to reduce 
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the possibility of reflections from the nozzle inflow. 

4. Inflow for the nozzle - characteristic conditions. Specifically, the Euler equations 
are linearized about the ambient state, assumed to hold far upstream in the nozzle, 
and the three incoming characteristics are imposed. 

p + pcu, v, Coop -p/c^ 

to be the values that they would have far upstream. It has been shown that this 
boundary condition is valid for the lowest propagating mode in the nozzle [9]. 

5. All other boundaries in the problem are artificial. Non-reflecting (radiation) 
boundary conditions are imposed to prevent spurious reflections from propagating 
into the interior. These boundary conditions are based on a far field expansion of 
the solution [22, 23]. 

A finite difference scheme is employed which is fourth order accurate in space and second 
order in time. The scheme is a generalization of the second order MacCormack scheme to 
allow higher order accuracy in space [24], The scheme is discussed in detail in [9, 15]. 

3 Results 

The results are presented in two major sections; (i) simulation of a jet without a nearby 
structure (i.e. free jet) and (ii) the jet with the nearby structure. In each case, the jet 
exit Mach number is Mj = 0.9 and the surrounding fluid is at rest. The nozzle utilizes a 
converging geometry with a contraction ratio of 2:1 over a longitudinal distance of four exit 
widths. Figure 1 shows a sketch of the computational domain for the case which includes 
the structure. 

The grids used in this study utilize stretching to enhance resolution in y of the shear layer 
and resolution in x near the nozzle exit. The t/-grid is mapped to conform to the contracting 
nozzle, while the £-grid remains Cartesian. This deviates slightly from orthogonality in the 
vicinity of the nozzle contraction but returns to a standard stretched Cartesian grid in all 
other regions. The computations have been extensively validated by grid refinement [1, 2] 
and, equally important, there are no spurious boundary reflections visible for the time interval 
considered in this paper (tc^lD < 100). 
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3.1 Free Jet Results 


For this case, non-reflecting boundary conditions are imposed on the upper boundary 
and the jet is centered in y. The origin of the x-coordinate is at the nozzle exit. The domain 
extends 45 D downstream, 30Z) upstream and 30D in both directions normal to the jet axis. 
The computed solution is symmetric around the jet axis to a close approximation, although 
symmetry is not imposed in the computation. A 581 x901 grid is employed. 

The structure of the pressure and vorticity field near the nozzle exit are considered first. 
Shown in figure 2 are contours of the instantaneous vorticity and pressure fluctuation, 

p(x, y , t) = p(x, y, t) - p(x, y), 

where p denotes the mean pressure, for a small region near the nozzle exit at tc^j D = 100. 
The sources simulated by the model are vorticity disturbances shed from the nozzle lip due 
to the near discontinuity in the velocity profile at the lip. The shedding is nearly periodic 
and the vortices convect primarily downstream. As they propagate, pressure disturbances 
associated with the vortices travel both upstream and downstream, ultimately leading to 
far field sound. The effect of the flow field inside the nozzle on the wavelength of the 
upstream propagating disturbances can be seen from the shorter wavelengths due to the 
slower upstream wave speed in the contracted section. 

Figure 3 shows the time history, spectrum and phase plot of the fluctuating pressure, p, 
at a point inside the jet flow field (x = 10D, y — y ax is = 0.8D). The PSD is normalized to 
0 db based on the maximum of the spectra. The figure shows that the spectral structure 
of the flow field p is composed of a fundamental frequency, /*, such that the Strouhal 
number, Sf = f*c ao j D ~ 0.25. The spectrum also shows two harmonics and one quarter 
subharmonic. This fundamental Strouhal number is close to observations of peak jet output 
(e.g., [25]). The precise value of 5 f * is certainly a property of the particular model used in 
the computations. However, this scaling for the shedding frequency has been confirmed for 
the model by computations with different jet Mach numbers in a limited range around 0.9 
and with different jet widths. Examination of p over an earlier time window (not shown) 
indicates a more broadband structure due to sound generated from the leading instability 
wave generated from the starter pulse. The predominant long time feature of the excited 
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jet predicted by the model is the nearly periodic shedding of vorticity from the nozzle lip. 
Computations with lower speed jets (not shown) indicate that this behavior becomes much 
more pronounced as the jet approaches a sonic exit velocity. 

The unsteady flow field in the jet generates acoustic disturbances which propagate into 
the far field as sound. Figure 4 shows a color representation of p over the entire computational 
domain at tCoo/D= 100. Although the flow is initially excited by a starter pulse, at the time 
for the figure this initial disturbance has propagated through and exited the computational 
domain. The remaining disturbances are generated due to instabilities of the jet shear layer. 
These disturbances convect downstream and spread with the jet. Analogous structures have 
been previously observed and computed in axi symmetric models [7, 9]. The figure also shows 
a series of waves propagating upstream into the far field together with a similar series of waves 
propagating downstream into the far field. The waves appear to be nearly periodic and are 
of similar spatial extent. Furthermore, they have the same frequency and wavelength. These 
waves are due to the nearly periodic shedding of vorticity from the nozzle lip (see figure 2) 
and the subsequent generation of pressure disturbances both inside and outside the jet. It 
can be seen from figure 4 that the upstream propagating waves appear to emanate from a 
point closer to the nozzle exit than the downstream propagating waves and also appear to be 
out of phase with the downstream propagating waves. This results in a region of interference 
between these two wave patterns along a ray directed approximately 100° from the jet axis. 

There is also a wave pattern within the contracting nozzle. The fluid field within the 
nozzle is nearly constant in y, leading to a pattern of nearly plane waves. The boundary con- 
dition inside the nozzle is rigid up to x = — 10 D and then is gradually changed to an absorbing 
boundary condition to help eliminate reflections from the inflow boundary. Animation of 
this data indicates that only upstream propagating waves are present in the nozzle. 

The far field pressure fluctuation is plotted in figure 5 as a function of fc^/D at three 
far field angles (0 = 30°, 90° and 120°) taken on a circle of radius 20 D. At all three points 
selected, p oscillates with a nearly constant period with a gradually varying amplitude. The 
period is close to that of the vorticity generation from the nozzle lip. This phenomenon is 
transmitted to the far field as sound. 

In addition to the time history, the PSD of p and phase plots are shown in figure 5. 
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The PSD is normalized to 0 db based on the maximum of the spectra over the three angles. 
The figure shows that the far field spectral structure is similar to the jet flow field spectra 
(figure 3). The fundamental frequency corresponds to ~ 0.25 and the two harmonics and 
one quarter subharmonic are exhibited as in the jet flow. The amplitude of the first and 
second harmonics is increased, relative to the fundamental, if the amplitude of the starter 
pulse is increased, confirming that these harmonics are generated via nonlinear interactions. 
Furthermore, the presence of low level subharmonics in the pressure spectra, again indicative 
of weak nonlinearities, are seen. Finally, it is seen that the spectra reach a maximum at 
9 ~ 30°, consistent with experiments and previous computations [2, 8, 9]. The relative 
harmonic content of the time signal increases as 9 increases. Thus there is a preferred 
upstream directivity for high frequencies (i.e. the harmonics). 

3.2 Jet /Structure Interaction 

Referring again to Figure 1, a wall, assumed infinite in extent and parallel to the nozzle, 
is located 3.5 D above the jet. The wall is assumed rigid, except for four flexible, aluminum, 
aircraft-type panels with clamped boundaries. Each of the panels is of length 3.5 D and 
thickness 0.01D. The panels are centered at xjD — 0, 3.75, 7.50, and 11.25, respectively. 
These are referred to as panels 1, 2, 3 and 4. Characteristic parameters of the panels are 
typical of aluminum. 

The grid in the jet domain is comparable to that described above for the free jet. In 
order to allow a greater distance between the jet and the artificial lower boundary the y- 
domain extends down 45Z). The panels are at the origin in the vertical direction. The 
Euler equations are also solved in the radiation domain. Using the fully coupled model, 
the unsteady disturbances in the jet are computed, which serve to excite the panels, the 
responses of the panels, and the radiation from the panels into both the jet domain and the 
radiation domain. 

A perspective plot of the unsteady flow and acoustic field can be seen in figure 6. In the 
jet domain the large scale structure propagating downstream is observed, trailed by other 
pressure disturbances in the jet. The generation of disturbances from the nozzle lip is clearly 
visible and occurs on a continual basis as time evolves. There is a series of ripples corre- 
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sponding to acoustic disturbances which radiate into the far field. The upstream radiation 
again appears to emanate from sources close to the nozzle exit, while the downstream radia- 
tion emanates from sources further downstream resulting in an interference pattern around 
0= 100°. The waves radiated downstream have higher amplitudes and the peak radiation 
occurs at around 30°. 

A series of outwardly propagating ripples, representing acoustic radiation from the panels, 
can be seen in the radiation domain. Since these ripples emanate from all four panels it is 
not possible to ascribe a precise geometric origin to them. Thus, instead of angles, far field 
locations in the radiation domain are designated by x/D values along a line at y/D = 25. 
The wavelength is comparable to the wavelength of the downstream radiated sound in the 
jet, i.e. the wavelength and frequency is controlled by Sf. The wave pattern in the radiation 
domain is more sharply defined that that observed in [1] for a lower Mach number jet. The 
results in [1] concentrated on earlier times and slower jets than the present computation. 
The largest overall level of sound radiation occurs at roughly 90° and the level is larger in 
downstream directions than in upstream directions. 

Features of the fully coupled flow/acoustics/ structural flow field are described next. The 
jet domain, panel response and radiation domain are examined. 

3.2.1 Jet Flow Domain 

In figure 7 the far field directivity is shown for both the free jet and the jet in the presence 
of a nearby wall. The data is given along an arc 20 D from the jet exit and for angles below 
the centerline. Far field pressures with the wall are similar to those described previously for 
the jet without any wall. Over the time interval considered, the differences are small. 

In figure 8, the incident pressure at the center of the panels is shown in both the time 
and frequency domains. Phase plots are also shown. The overall level of p increases with 
the downstream distance of the panels. The level of p on panel 1 (partially upstream of 
the nozzle exit) is reduced nearly an order of magnitude from the loading on panel 4. Since 
panel 4 is at a lower angle to the jet axis, this is in agreement with the directivity results 
in figure 7. The frequency spectra of the incident pressure is dominated by the fundamental 
frequency, and its harmonics. The relative harmonic content of the incident p for panel 1 
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is significantly increased as compared to the other (more downstream) panels, due to the 
preferred upstream beaming of high frequencies. 

3.2.2 Panel Response 

The vibrating velocity, u, at the center of each panel is plotted in figure 9, together 
with the corresponding PSDs and phase plots. The panels were held fixed until the leading 
pulse had passed. This eliminates any effect of the pulse in the panel response, but does 
lead to an impulsive start to the panel motion. For this reason, only the long-time response 
of the panels is investigated. The panel response increases with distance until reaching its 
maximum on panel 3 and then decreases slowly with distance. The relative subharmonic level 
is increased as compared to the incident pressure. This may be due to enhanced receptivity 
of the panels to low frequency excitation. The predominant response of the panels, however, 
is at the jet fundamental frequency, S*, and its harmonics. The relatively large harmonic 
contribution to the response of panel 1 is indicative of the preferred upstream beaming of 
high frequencies, although the overall response is smaller than that of the more downstream 
panels. 

The observation made in [1, 2] that the panels act as narrow filters, converting a broad 
band input into a narrow band response, may no longer be relevant for jets dominated by 
the long-time vortex shedding considered here. In this case the panels respond clearly at the 
frequencies of the loading, somewhat enhancing the subharmonic response. 

3.2.3 Radiation Flow Domain 

The panel radiation is considered next. The pressure, p, is computed along the line 
y = 25 D in the radiation domain. Shown in figure 10 are time histories, PSDs and phase 
plots for four different ^-locations along this line. It can be seen from the figure that the 
overall level of radiation is reduced as x decreases, i.e. in the upstream direction in terms of 
the jet flow field. The PSDs indicate that for downstream locations, there is peak radiation 
at S*, the first two harmonics and the subharmonic. For upstream locations, a low level, 
relatively broad band spectrum is found. Thus the radiated sound is beamed primarily in 
vertical and downstream directions and is controlled by S'*, i.e. the vortex shedding frequency 
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of the jet. 

In figure 11 the overall intensity is computed as a function of x. This figure illustrates the 
strong vertical and downstream beaming of the panel radiation. It also illustrates a zone of 
relative silence at about xjD — 15 where the intensity is reduced by about 15 db from nearby 
locations. This zone can also be seen in the flow field in figure 6 and is due to destructive 
interference between radiation from the various panels. No such zone was observed at lower 
jet Mach numbers [1, 2], where the resulting radiation was more heavily influenced by the 
starter pulse and the leading instability wave generated sound in the jet. 

In a real jet, instabilities can be generated by bursts within the jet, generating large scale 
structures, as simulated by the starter pulse and the intermediate time behavior of the jet. 
Also, instabilities can be generated in a more continual fashion via vortex shedding, as is 
done here, by considering the long time evolution of the model. The model allows for a 
separation of these effects so that they can be studied separately. 

4 Conclusions 

A numerical simulation has been conducted of a nearly transonic jet exhausting into a 
fluid at rest both with and without the presence of a nearby structure. The structure is rigid 
except for four flexible aircraft-type panels. These panels respond to acoustic forcing from 
the jet and radiate sound into both the jet domain and a domain on the other side of the 
panels. The simulation utilized a full coupling between the Euler equations (solved on either 
side of the panels) and the nonlinear beam equation resulting in a full coupling between the 
flow, acoustics and panel response and radiation. The results are in qualitative agreement 
with experimentally observed phenomena and with previous computations indicating the 
present model is capable of capturing the important features of jet flow/acoustic/structural 
interactions at nearly sonic jet velocities. 

In particular, the results indicate the following properties of the long time response of 
excited, high subsonic Mach number jets: 

1. There is a nearly periodic shedding of vorticity from the nozzle lip. This shedding 
is associated with pressure disturbances which lead to far field sound. 
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2. The shedding frequency scales with Strouhal number. S* is near 0.25 which is 
near the peak of many observed far field radiation patterns. 

3. The relative harmonic content of p increases with upstream angle. 

4. The far field sound peaks near 30°. 

All of the above are consistent with the features of previously observed real jets. Furthermore, 
the panel response and radiation results indicate that: 

1. Loading and response of the panels generally increases with downstream distance. 

2. The frequency spectrum of the panel responses is dominated by the jet shedding 
frequency, S*, its harmonics and subharmonics. 

3. The filtering action of the panels previously reported [1, 2] is less evident in the 
present work, where the panel forcing is dominated more by the periodic shedding 
from the nozzle lip and is already of a narrow band nature. 

4. The relative level of harmonics in the panel loading and response is increased for 
upstream panels. 

5. Radiation is primarily in vertical and downstream directions. 

6. The upstream radiation is small compared to downstream radiation and does not 
appear to reflect the jet shedding frequency. 

7. Zones of relative silence occur due to destructive interference. 

Finally, the code developed is computationally efficient, running at over 190 MFLOPS on a 
Cray-YMP. This produces a complete solution of the fully coupled flow/acoustic/structural 
interactions in approximately 6 CPU hours. 
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Figure 1: Computational Domain 
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Figure 2: Near field vorticity and pressure contours 
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Figure 3: Jet flow field time history, spectrum and phase plot at x/D = 10 and 
{y ~ V axis) l -D — 0.8 
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Figure 4: Contours of p for free jet 
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Figure 9: Panel velocity time histories, spectra and phase plots 
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Figure 10: Radiation domain far field time histories, spectra and phase plots at 
x/D = —5, 10,25, and 40 along the line y/D = 25 
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Figure 1 1 : Radiation domain directivity of overall sound pressure level along the line y = 25D. 
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